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Abstract 



The flavor conversion of supernova neutrinos can lead to observable signa- 
tures related to the unknown neutrino parameters. As one of the determinants 
in dictating the efficiency of resonant flavor conversion, the local density pro- 
file near the MSW resonance in a supernova environment is, however, not so 
well understood. In this analysis, variable power-law functions are adopted 
to represent the independent local density profiles near the locations of reso- 
nance. It is shown that the uncertain matter density profile in a supernova, 
the possible neutrino mass hierarchies, and the undetermined 1-3 mixing angle 
would result in six distinct scenarios in terms of the survival probabilities of v e 
and v e . The feasibility of probing the undetermined neutrino mass hierarchy 
and the 1-3 mixing angle with the supernova neutrinos is then examined using 
several proposed experimental observables. Given the incomplete knowledge 
of the supernova matter profile, the analysis is further expanded to incorpo- 
rate the Earth matter effect. The possible impact due to the choice of models, 
which differ in the average energy and in the luminosity of neutrinos, is also 
addressed in the analysis. 
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I. INTRODUCTION 



Definite evidences from recent experiments [1-7] have confirmed that the neutrinos are 
massive and mixed. Part of the intrinsic neutrino properties, such as the neutrino mass 
hierarchy and the 1-3 mixing angle, however, remain undetermined. The need for a better 
knowledge of the neutrino physics has motivated plenty of phenomenological analyses on 
whether and how the unknown neutrino properties can be reliably probed with the next 
generation experiments involving neutrinos from both the astrophysical and the terrestrial 
sources [8-14]. 

The core-collapse supernovas represent a unique type of neutrino source in that they emit 
neutrinos of all three flavors with a characteristic energy range in a time scale distinct from 
that of the neutrinos emitted from the sun, the atmosphere, and the terrestrial sources. The 
rich physical content involved in the environment further makes the supernova an important 
platform for the study of neutrino properties. In particular, the neutrino burst from a su- 
pernova has long been considered as a promising tool for probing the undetermined neutrino 
parameters, despite the fact that the complexity arising from the unavoidable astrophysical 
uncertainties may result in ambiguity in the interpretation of the observed events. 

The standard MSW formulation [15,16] suggests that as the supernova neutrinos prop- 
agate outward, the variation of density profile can have an impact on the resonant flavor 
conversion. One of the difficulties in probing the unknown neutrino properties with the 
supernova neutrinos arises from lack of knowledge to the matter density profile in a super- 
nova. The inverse power-law density, p ~ r~ 3 , of the progenitor star is usually adopted in 
the literature as the density distribution through which the supernova neutrinos propagate 
outwards. However, whether the global density distribution p ~ r n with a fixed power can 
provide a justifiable connection between the dynamics of resonant flavor conversion and the 
expected neutrino events at the detectors remains uncertain. In addition, one may easily 
question the validity of this over-simplified model if, among others, the time-dependent mat- 
ter distribution driven by the shock waves and the associated effects (see, e.g., Ref. [17-22]) 
are considered. It is shown [23] that near the MSW resonance, even a relatively small de- 
viation of the power from n = — 3 could have a sizable impact on the interpretation of the 
neutrino parameters. Furthermore, no evidence suggests that the density distributions near 
both the higher resonance and the lower resonance should be described by a similar density 
profile. 

Another major source of astrophysical uncertainty in analyzing the supernova neutrinos 
would be the spectral parameters. Distinct predictions to the average energy and the lu- 
minosity of neutrinos exist among models. Thus, the possible impact due to the uncertain 
spectral parameters also needs to be addressed in an analysis. 

Independent and variable power-law functions are assumed in this work as the matter 
profiles near the locations of MSW resonance. In addition, physical observables derived from 
the expected neutrino events at a water Cherenkov detector are proposed. It is the aim of 
this work to analyze the outcomes due to the uncertain density profile and the choice of 
spectral models, and to investigate whether and how the experimental observables would 
provide hints for the determination of the unknown neutrino parameters. The analysis is also 
extended to include the Earth matter effect. This work is organized as follows. Section II 
summarizes the general features of the neutrino fluxes emitted by a core-collapse supernova 
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and that of the relevant detection processes at a water Cherenkov detector. In section III, the 
variable power-law function is proposed to account for the possible uncertainty in the matter 
profile near the MSW resonance. In addition, observables are constructed and examined. 
Focus has been placed on whether and how these observables would resolve the various 
scenarios arising from the uncertainty in the matter profile, the possible neutrino mass 
hierarchies, and the undetermined mixing angle $13. In Section IV, the general formulation 
is expanded to incorporate the regeneration effect due to the Earth matter. In Section V, the 
observables proposed in Section III are re-analyzed with the Earth matter effect included. 
Section VI is devoted to analyzing the impact due to the choice of spectral models. We then 
summarized this work in Section VI. 



II. SUPERNOVA NEUTRINOS AT THE DETECTOR 

A core-collapse supernova emits neutrinos of all three flavors on a typical time scale of 
roughly ten seconds. Within a factor of two, the total luminosity is approximately equi- 
partitioned into each flavor. The mean energy of each neutrino flavor is determined by its 
interacting strength with matter. It is believed that the neutrino spectra are not exactly 
thermal; they are usually modeled by the pinched Fermi-Dirac distributions. The flux of 
flavor v a reaching the Earth surface can be written as 

r p2 

a 47rr 2 T^F 3 { Va ) d( E ^ T ^'^ + 1' K) 

where L a is the luminosity of the neutrino flavor u a , r is the distance to a supernova, T a is 
the effective temperature of u a inside the respective neutrinosphere, E a is the energy of u a , 
and rj a is the pinching parameter. The normalization factor F 3 (i] a ) is given by 

F 3 ( Va ) = -6Lu(-e^), (2) 

with the polylogarithm function Li n {z) defined as 

^»(*) = f;£ (3) 
k=i K 

Note that the temperature T a is related to the mean energy (E a ) by 

{Ea) ~ 3 L* 3 (-^) a ' { ) 

As an illustration, the spectral parameters published by the Lawrence Livermore (LL) 
group [24] will be adopted in the analysis: (E Ue ) : (E 9e ) : (E Ux ) ^4:5:6 with (E Pe ) = 15 
MeV, and L v jL Vx = 2.0, L D jL Vx = 1.6, where v x = u^, v T , z/ M , and v T . The pinching 
parameters i] Ue = i] Pe = r\ Vx = 3.0 are assumed. As for the neutrino parameters, we choose 
8m 2 21 = 7.0 x 10" 5 eV 2 , |5m^| = 3.0 x 10~ 3 eV 2 , and sin 2 29 12 = 0.81. Note that the ongoing 
work by the Garching group [25,26] predicts distinct neutrino fluxes and average energies 
from that of the LL group. The comparison will be performed and discussed in Section VI. 
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In a water Cherenkov detector, such as the Super-Kamiokonde (see, e.g., [1]), the up- 
graded Hyper-Kamiokonde [27], and other facilities that are under planning or construc- 
tion [28], the event rate induced by the isotropical inverse /3-decay process 

V e + V — > n + e+ (5) 

dominates over that of other channels. In this analysis, the isotropical CC processes: 

u e + O^F + e~, (6) 



v e + O -> N + e+, (7) 

and the highly directional events induced by u a + e" (a = e, p,, r) and V a + e" scattering 
will also be included. 



III. CONSTRUCTING OBSERVABLES FROM THE NEUTRINO SIGNALS 

A. The formulation 

The neutrinos would likely go through multiple flavor transitions in a supernova at 
distinct density scales. A single power-law function, p ~ r~ 3 or p ~ r n in general, is 
usually adopted in the literature as the density profile of the supernova matter. When the 
neutrinos travel through a complex environment such as the matter of a supernova, however, 
it would seem more plausible in an analysis to allow independent and variable local density 
profiles near the resonant locations. 

If the matter distribution is modeled as layers of matter along the path of neutrinos in 
a specific direction, a simple power-law function can be adopted as the local profile within 
a layer: 

p(r) k = c k r n \ (8) 

where k stands for the /c-th layer, r nk describes the matter profile within the k-th layer, and 
Cfc denotes the magnitude of the component r Uk . The density profile near the location of 
resonance can then be written as 

p(r) fco =c fco r"*o, (9) 

where the fc -th layer stands for the resonance layer of the flavor conversion. Note that in 
this approach the power can vary from layer to layer to account for the uncertainty of profile 
originating from various types of sources. The thickness of the layers, which do not take part 
in the formulation of this analysis, can be uneven and unrestricted, with the assumption that 
a specific power n properly describes the density profile within a specific layer. Furthermore, 
the variable power and thickness also allow the deviation of the matter density from spherical 
symmetry and the inhomogeneity of matter distribution along different directions. 

One may also employ functions of more than one polynomial term. The extra polynomial 
terms would provide more details about the density profile in a global scale. However, the 
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trade-off is that they also introduce more undetermined coefficients for the extra terms. 
Since one only concerns the local density profile near the MSW layer, the detailed shape 
and variation of the density in a global scale would be unimportant, and the extra terms with 
their undetermined coefficients would become redundant in the analysis. Thus, the power- 
law function with an adjustable n would be suitable and much simpler for the analysis on 
how the variation of local density profile, rather that the global profile, would impact the 
MSW resonance. 

In the following, one may simply focus on the density shapes of matter in the resonance 
layers: p(r) = C\r ni for the lower resonance, and p(r) = Chr Uh for the higher resonance. The 
variation of n/ and n h would be expected to alter the efficiency of neutrino flavor transition. 

Along this approach, the crossing probabilities at the resonance can be modified as 



exp(-f lk F k ) - expj-f 7fcgg%) 



(10) 



where k = h or I, 9^ is the mixing angle, F k is the correction factor to a non-linear profile, 
and the adiabaticity parameter 7^ takes the form 

_ 1 M% sin 2 2^ cos2% 1 

lk ~ 2|n fe | ( E v ] ( cos2% K 2 v / 2G F ^c fc ) ' ( } 

where Srrifj = \mf — m 2 \, Gf is the Fermi constant, E v is the neutrino energy, Y e is the 
electron number per baryon, m n is the baryon mass, and the scale c k varies very weakly 
with r over the range 10 12 g/cm 3 < p < I0~ 5 g/cm 3 . Note that P h = P h , and that P t can be 
obtained directly from Pi by replacing % with 7r/2 — 

The probabilities that v e and u e survive the flavor conversion are given respectively by 

P nor ~ U 2 el PtP h + U 2 e2 (l - P{)P h + UUl - P h ), (12) 

P^-UKl-P^ + U^Pu (13) 

for the normal hierarchy, and 

P inv ^^\-P l ) + Ul 1 P u (14) 

P inv ~ XJl 2 P { P h + U 2 el (l - P^Pn + U 2 3 (l - P h ), (15) 

for the inverted hierarchy. Eqs. (10) and (11) imply that the variation of neutrino energy 
plays an insignificant role [23], as compared to the variations of n k and the 1-3 mixing angle, 
in affecting the survival probabilities. In addition, the two extreme values of Ph. Ph ~ 1 
(non-adiabatic) and Ph ~ (adiabatic), are given by the conditions: Ph ~ 1 if g < 1, and 
P h ~ if g > 1, where the function g relates n h and 9 13 : 

9 = 9(n h ,9 13 ) ~ 2 ' 3 , X ^ (37.6 x lO" 30 x cos2^ 13 )^^^. (16) 
\n h \ cos26'i3 
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Furthermore, Eqs. (10) and (11) also suggest that the arbitrary n\ gives rise to only two 
distinct values of P\ (Pi): Pi ~ (Pi ~ 0) for rii > —5, and Pi ~ cos 2 9 12 (Pi ~ sin 2 # 12 ) for 
n t < -5. 

As summarized in Table 1, the variation of n k (k = /, h), the undetermined #i 3 , and the 
possible mass hierarchies would lead to six distinct scenarios in terms of the survival prob- 
abilities. In the first row of Table 1, the quantities in square brackets denote the variables, 
i.e., Pi and Pi are determined by the variable ni; Ph is determined by nu and g; and the total 
survival probabilities, P and P, are determined by all rii, rih, g, and the mass hierarchy. The 
last column labels the six different scenarios in terms of the possible combinations formed by 
P and P. Note that the variation of n k would impact the probabilities differently according 
to whether n k > —5 or n k < —5. Furthermore, the outcomes due to p ~ r~ 3 , which is 
usually assumed in the literature, correspond to those due to (rii > —5, rih > —5) in Table 
1. 

Some qualitative observations are outlined as follows, (i) No information on the mass 
hierarchy or on 9 n is available if either scenario I or scenario IV is observed, (ii) If either II 
or III is observed, then the lower bound on #13 would be available, (iii) Likewise, if either V 
or VI is observed, then the lower bound on #13 would be available. 

The neutrino events may be categorized into isotropical ( Eqs. (5), (6), (7)) and di- 
rectional (neutrino scattering) events at the water Cherenkov detector. The expected time- 
integrated number spectra for both the isotropical and the directional events are shown in 
Figs. 1(a) and 1(b), respectively. It can be seen from Fig. 1(a) that three groups of nearly 
degenerate spectra would form for the isotropical events: (I,V), (II, IV), and (III, VI). Thus, 
the isotropical events alone fail to provide all the needed clues for separating the six scenar- 
ios. On the other hand, while the directional event rates shown in Fig. 1(b) amount to only 
a tiny fraction of the isotropical ones, they could provide the rest of the clues for lifting the 
above degeneracy. As a result, combining both the isotropical and the directional events is 
inevitable in the construction of feasible observables. 

The v e and v e fluxes arriving at the surface of Earth, F e and F s , are related to the original 
ones, F® and F® , through the survival probabilities P and P for v e and v e , respectively: 

F e = F e °+(l-P)(F°-F°), (17) 

F g = F°+(l-P)(F°-F°). (18) 

The expected total event number N a induced by the flavor u a at the detector is given by 

N t 

Na = 4^J T,l F ^¥E, (19) 

where N t is the number of targets at the detector, r is the distance to the supernova, and 
c"qi = (?ai(E) is the cross section for v a in a particular reaction channel i. The detection 
efficiency is assumed to be one. Note that the regeneration effect due to the Earth matter 
will be expected to alter the survival probabilities, and thus the observed event rates. We 
shall first concentrate on building the observables without considering the Earth regeneration 
effect. A more detailed analysis incorporating the Earth effect will follow. 
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B. Observables derived from neutrino signals 



Statistically, the predominant isotropical channels can provide the most satisfying events 
for the analysis. Fig. 1(a) suggests that one may obtain the optimized results by examining 
R 1 , which is defined as the ratio of the event rates in the medium and low energy range 
(E < 40 MeV) to that in the high energy tail (E > 40 MeV): 

N U E<A0MeV) 

N iso (E > AOMeV) 1 ; 

This ratio can be reduced to 

. ..1 + 0.30P - 0.008P 

R 1 ~ 2.62 = . 21 

1 - 0.88P - 0.028P v ; 

This expression leads to the following possible values for R(, where J labels the six scenarios: 
(i) R[ = R\ ~ 6.6. (ii) P{ 7 = R[ v ~ 8.8, (iii) R[ H = R\ ! ~ 2.6. It is seen that the values 
of Pi for the group (III,VI) are smaller than the others by a factor of 2 ~ 3 and can be 
easily singled out. 

From a different point of view, one notes that the inverted mass hierarchy implies the 
following possible results: (i) Pi ~ 6.6, which predicts n/ < —5, rih < —5, while the 
information of 13 is unavailable, (ii) Pi ~ 8.8, which predicts rii > —5, rih > —5, g < 1, 
and thus the upper bound sin 2 9 13 ~ 10~ 2 . (iii) Pi ~ 2.6, which predicts ni > —5, g > 1, and 
thus the lower bound sin 2 13 ~ 10~ 2 . On the other hand, the following would be possible 
if the mass hierarchy is normal: (i) Pi ~ 8.8, which predicts n\ > —5, while rih and #13 are 
undetermined, (ii) Pi ~ 2.6, which only predicts rii < —5. 

To establish the P-sensitive observables that resolve the degenerate scenarios in the same 
group, we now examine the z/ e -dominating directional events. To optimize the results, one 
may define the ratio of the directional events for E < 25 MeV to that for E > 25 MeV as 
P 2 : 

N dtr (E < 25 MeV) 
2 N dir (E > 25MeV)' 1 ' 

This ratio can be reduced to 

„™ 1 + 0.23P + 1.19P 
R - ^ 98 l-0.08 6 P-0.3 6 f (23) 

which leads to the following: R{ ~ 2.46, B% ~ 1.22, R\ u ~ 1.45, R 1 / ~ 1.76, R\ ~ 1.17, 
and R 1 / ~ 2.13. 

It is seen that the six scenarios lead to distinct observable results in terms of the com- 
bination (Pi, P 2 ). The observable Pi would reduce the six scenarios into three degenerate 
groups: (I,V), (II, IV), and (III, VI). Observing P 2 would further remove the degeneracy be- 
tween the scenarios in each group. Thus, each scenario can be characterized by a unique set 
of observables (Pi, P 2 ), which should act as a feasible reference in removing the degeneracy 
among the six scenarios. 

In searching for useful observables other than Pi and P 2 , one notes from Fig. 1 that 
in terms of the height and the width, the possible spectra exhibit the most pronounced 
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differences in the region near the peaks. The ratio formed by the peak value of a specific 
spectrum (F max ) and the width at half peak value (T^) deserves some attention. One may 
introduce the double ratio involving both the isotropical and the directional events: 

j£ _ (Fmax /T 1/2) iso (24) 
(Fmax /r 1/2) dir 

The six scenarios lead to the following possible values: R{ ~ 16.7, ~ 49.5, R{ n ~ 22.6, 
R 1 / ~ 35.5, R% ~ 40.4, and R^ 1 ~ 11.9 

Furthermore, one may define another observable P4 involving both the isotropical the 
directional events in the medium energy range, 13MeV < E < 27MeV: 

N ios (13MeV < E < 27MeV) 
4 ~ N dir (13MeV <E< 27MeV) ' 1 ' 

which leads to 

1 + 1.05P + 0.0013P 
R ' = 19 86 1+0.18P + 0.59P ' P6) 

and R{ ~ 22.3, R 1 / ~ 31.3, Pf J ~ 17.0, Pf y ~ 27.2, P y ~ 29.5, and R\ T ~ 14.7 

If the neutrino energy spectrum is reconstructed, one may examine yet another observable 
i? 5 , which involves the average energies of both the isotropical and the directional spectra: 
R 5 = (E) iso /{E) dir . The combined average energy (E)i is defined as 

{ )l = — y^T ' (27) 

where i stands for isotropical or directional, is the time-integrated number spectrum, 
N a i is the event number, and the summation over a includes all the possible isotropical or 
directional channels induced by v a . Further reduction of R^ yields 

(1 - Q.37P - 0.018P)(1 + 0.71P + 0.46P) 
5 ' (1 -0.023P-0.014P)(1 - 0.002P-0.043P)' 1 J 

and the values of P5 for the six scenarios are given by: R\ ~ 1.27, R^ 1 ~ 0.99, R^ 11 ~ 1.41, 
R J 5 V ~ 1.10, P^ ~ 1.04, and P^ 7 ~ 1.56. 

The uncertainty arising from the density profile can leave observable signatures on the 
neutrino events through the resonant flavor transition. For a given set of input spectral pa- 
rameters, the proposed experimental observables may be combined to form useful references 
for removing the degeneracy among scenarios. However, it should be noted that whether the 
unknown neutrino parameters could be definitely determined remains as an open question 
as other possible sources of uncertainty may blur or even wash out the signatures. We shall 
discuss this issue in sections VI and VII. 



IV. INCLUDING THE EARTH MATTER EFFECT 

The neutrinos emitted from a supernova arrive at the Earth surface in mass eigenstates. 
Since the matter density of Earth is approximately the same order of magnitude as that 
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near the lower resonance in a typical supernova, the neutrinos could oscillate again if they 
cross Earth prior to the detection. This regeneration effect in the Earth matter would be 
expected to alter the neutrino fluxes arriving at the detector. The observable effects due to 
this regeneration process have been intensively studied [29]. It is realized that in general 
the Earth effects are signaled by an oscillatory modulation on the neutrino spectrum that is 
expected when the regeneration effect is missing. Given the uncertainty in matter density 
profile of a supernova, we intend here to further examine the feasibility of the proposed 
observables under the influences of the Earth regeneration effect. 

The regeneration effect depends on several factors: the neutrino parameters, the neutrino 
energy, the density shape of Earth matter, and the location of the detectors (or, the incident 
direction of the neutrinos). The v e and v e fluxes arriving at a detector D can be expressed 
as P e D and PP, respectively: 

F e D = F e °[P,(P, + P 2e (l - 2P t ))] + P°[l - P h (P t + P 2e (l - 2P l ))], (29) 



F? = P°[l - (P t + P 2e (l - 2P))] + P°[P + P 2e (l - 2P)] (30) 
for the normal hierarchy, and 

F e D = P e °[P + P 2e (f - 2P)] + P°[l - (P + P 2e (l - 2P))], (31) 



F? = F°[P h (l - (P + P 2e (l - 2P))) + P°[l - P h (l - (P + P 2e (l - 2P l )))\ (32) 

for the inverted hierarchy, where P 2e (P 2e ) represents the probability that a z/ 2 (z/ 2 ) arriving 
at the Earth surface is detected as a v e (z/ e ) at the detector. 

The Earth matter density encountered by the electron neutrinos is given by V(x) = 
V / 2Gi?iV e (a;), where Gf is the Fermi constant and iV e (:r) = p{x)/{m p + m n ), with p(x) the 
density profile, m p the proton mass, and m n the neutron mass. The probability P 2e is given 
by [30,31]: 

P 2e = sin 2 6 12 + \ sin 2 2^ 12 H dx[V(x) smcf) x ^ Xf }, (33) 

where Xq is the entry point, Xf is the exit point, and x is the distance into the Earth. The 
function (p a ^b is defined as 

b A(x)dx, (34) 



where 



A(x) = ^|i v /[ cos2 ^-e(x)] 2 + sin 2 2^ 12 , (35) 



2E 

and the parameter e(x) is given by e(x) = 2EV \x) / '5m 21 . Note that 



e(x) » 7.575 x 10' ^ ' 2 ^f' (36) 
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and 

^« 2.543 x 10 3 ^f#^(l/M. (37) 
2E E/MeV y 1 ' y ' 

The density profile of the Earth matter can be approximated by a simple step func- 
tion [32,33]: p(x) 5.0 g/cm? for |i? e < R < i? e (mantle), and ~ 12.0 g/cm 3 for 
i? < ^i? e (core), where is the Earth radius. For a specific detector, the incident angle of 
the neutrino determines the path length through Earth and the magnitude of P 2e . Because 
of the chosen density model, it is necessary to calculate P 2e separately for both < ip < i\) a 
(mantle only) and ip a < ip < tt (mantle+core+mantle), where ip is the incident angle in 
general, and ip a denotes the incident angle when the path is tangent to the sphere of radius 
R = \R(q, at which radius the density changes. 

For < ip < ipa, as shown in Fig. 2(a), the probability P 2e takes the form 

P 2e = sin 2 # 12 + isin 2 2# 12 [e(x)(l - cos(LA(x)))], (38) 
where the total path length L inside the Earth is given by 



L = {d- Rq) cos^ + yjR% - (i? ffi - d) 2 sin 2 6 12 , (39) 

with d the depth of the detector. 

For ip a < tp < n, the path inside the Earth consists of three segments: x a b, Xb c , and x c d, 
as shown in Fig. 2(b): 



/ R 2 
x ab = L- (i? e - d) cos^ - J (Re - d) 2 cos 2 ^- [(i? e - d) 2 - — ], (40) 



R 2 

x bc = 2J(R (S - d) 2 cos 2 ^- [(i? e - d) 2 - — ], (41) 



R 2 

x cd = (d - R e ) cos^ - y (i? e - d) 2 cos 2 ^ - [(i? ffi - d) 2 - ~^\- (42) 
The probability P 2e in this case can be written as 

P 2e = sin 2 9 12 + X - sin 2 26 l2 x J, (43) 

where 

/ = ^(^[i _ cos(A m x a6 )] + e c (x)[l - cos(A c z 6c )] + e m (x)[l - cos(A m x cd )], (44) 

and the superscripts m and c stand for mantle and core, respectively. It is clear from the 
oscillatory terms that a longer path length in matter would give rise to a higher frequency 
of P 2e . Note that for the probability P 2e , the parameter e(x) = 2EV(x)/5ml 1 is replaced by 
-e(x): 

1 f x f f x f Ath? i 

P 2e = sin 2 9 12 - - sin 2 20 l2 / dxV (x) sin( / dx— J(cos20 12 + e(x)) 2 + sin 2 29 12 ). 

2 Jxq Jx 2E 

(45) 
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V. REEXAMINING THE OBSERVABLES WITH THE EARTH MATTER 

EFFECTS INCLUDED 



For one of the six scenarios J, the expected neutrino flux of flavor a (a — e, e) at 
the surface of Earth, F£, and at the detector, ,J , are determined by the type of mass 
hierarchy, the values of Ph, Pi (Pi), and the original fluxes F® for a given model. We list the 
expressions of F^ and F® yJ in Appendix A. The effective survival probabilities for v e and 
v e at the detector, P DyJ and P D ' J , can be read off from F® ,J and F®' J ', respectively: 

pDJ = cos 2 Q 12 + p 2e(1 _ 2 cos 2 9 12 ), 
P D ' n = sm 2 e 13 < 1, 

P D ' IH = P2e, 
P D > IV = P 2e , 

P D ' y = sin 2 13 <l, 
pD,vi = cos 2 9l2 + p 2e (! _ 2 cos 2 12 ), (46) 

and 

pDJ = cos 2 12 - P 2e (l - 2 sin 2 12 ), 
pA" = 1 - P 2e , 

pDJII = gin 2 < 1; 

pD,/y = ! _ p 2e? 

P^.V- = cog 2 _ p 2e(1 _ 2 sin 2 ^ 2 ) 5 

P D ' VI = sin 2 < 1. (47) 



In principle, the flux variation due to the Earth matter effect can be conveniently illus- 
trated by the ratios R e = (F e D — F e )/F e and i? g = (P g D — Fg)/Fe, for z/ e and z/ e , respectively. 
The ratios i?g and i?g, as given in Appendix B, suggest the following: 

(i) The values of R^ 1 , R\ , R 1 / 1 , and R^ 1 vanish, i.e., including the Earth matter effect 
would not alter the resultant v e flux for II and V, and the resultant v e flux for III and VI. 
This can be understood by the fact that the crossing probability at the higher resonance 
approaches zero (Ph ~ sin 2 #13 -C 1) for each of the above scenarios, as indicated by Table 
1. If the neutrinos travel through the Earth matter before reaching the detector, Table 1 
and Eq.(46) imply that P = P D)J = sin 2 9^ <^ 1 for scenarios II and V. Likewise, Table 1 
and Eq.(47) imply that P = P DyJ = sin 2 13 < 1 for scenarios III and VI. 

(ii) The two ratios, Rl 1 and R 1 / , are identical because of P D ' n = P D ' IV = 1 - P 2e . In 
addition, one may also conclude that Rl 11 = R™ , Rl = R V J , and R\ = R\. 

(iii) The scenarios in each degenerate pair, (II, IV) or (I, VI), will not be distinguished 
from one another by the observation of v e events alone even if the Earth matter effect is 
included. Thus, constructing the observables that are sensitive to the v e events would be 
beneficial to the analysis. Likewise, the z/ e -sensitive observables would be needed to lift the 
degeneracy between III and IV, as well as that between I and IV. 

In summary, the values of R e and R s may act as a useful qualitative and quantitative 
reference in analyzing the distinct scenarios. For example, scenario II and scenario IV can 
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be distinguished by examining whether the Earth matter effect is missing ( scenario II) or 
observed ( scenario IV) in the v e flux. Likewise, scenario III and scenario IV can be separated 
by examining whether the Earth matter effect is observed (scenario IV) or not (scenario III) 
in the v e flux. Similar qualitative observations can be applied to separating I and V, or I 
and VI. 

We now reexamine the expected ratio Pf = NP (E < 40MeV)/N° o (E > 40MeV) at 
a detector when the Earth matter effect is included. The impact on Pf from the Earth 
effect can be examined through the plot Pf = Ri(ip) in Fig. 3. As expected, scenarios III 
and VI can be easily singled out by their low values, Pf ~ 2.6. Note that the values of 
Pf for III and VI are nearly independent of the incident angle of neutrinos. In addition, 
a jump of Pf occurs at tp ~ n/2 for all the scenarios. This is realized by the fact that 
the length of neutrino path in mantle begins to increase significantly when ip > n/2. It is 
also seen that an even more significant jump occurs at ip ~ 5n/6, beyond which angle the 
neutrino flux starts to pass through the Earth matter that is formed by the combination 
mantle+core+mantle before reaching the detector. Note that the values of Pi for the four 
scenarios: I, II, IV, and V, become almost indistinguishable when if) > 57r/6. 

As discussed in Section III, the neutrino events from the P-sensitive channels are needed 
to remove the degeneracy that the observable Ri fails to resolve. Thus, the P-sensitive 
ratios P 2 = N^ ir (E < 2hMeV) / N^ ir (E > 25MeV) as functions of ip for all the six scenarios 
are reexamined and shown in Fig. 4. 

To explain the obvious step-up behavior of III and IV, and the step-down behavior of 
I, II, V, and VI in Fig. 4, one first notes from Eq. (43) that P2 e increases with I . While 
the path length determines the frequency of the oscillatory term, the mean value of P2 e is 
determined by the sum e m -|_ e c _|_ e m for 57r/6 < ip < ir, and by e m for < ip < 57r/6. Since the 
parameter e is proportional to the density, so that P 2e increases with the density. Likewise, 
Eq. (45) implies that p2 e decreases when the density increases. It then follows from Eq. (46) 
that P D)I and P D ' VI decrease when the density increases (1 — 2 cos 2 9\2 < 0), while P D > IU 
and P D > IV increase when the density increases. On the other hand, Eq. (47) implies that 
all P D,! , P D ' n , P D,IV , and P D < V increase when the density increases. Given the conditions 
that P D ' H , P Dy , P D > m , and P D ' VI are vanishing, and that the coefficients of the P-term 
dominate over that of the P-term in both the denominator and the numerator of Eq. (23), 
the stepping down behavior of Pf ' 7 , R®' 11 , R^ ,V \ and R^ Vl ' , as well as the stepping up of 
R^ IU and R? IV at greater become apparent. 

It should be pointed out, as mentioned earlier, that the Earth matter effect is manifested 
by an oscillatory modulation on top of the neutrino spectra that would have been observed 
without the regeneration effect. Practically, it would be very difficult, if not impossible, to 
extract from R® the information that is significantly different from that obtained by P3. 

The ratio of the isotropical events to the directional events, Rf = N® a (l3MeV < E < 
27MeV)/N^ r (13MeV < E < 27MeV), and the ratio formed by the average energy of the 
isotropical to that of the directional events, Pf = (E)f so /{E)® ir , are shown in Fig. 5 and 
Fig. 6, respectively. 

It is seen that the six scenarios are in principle distinguishable by examining Pf , R21 
Rf and R®. However, as can be seen from the figures, the distinct scenarios become less 
distinguishable at a larger incident angle, 5n/6 < ip < 71. Given the potential theoretical and 
experimental uncertainties, ambiguity may arise in interpreting the proposed observables if 
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the Earth matter becomes pronounced. 



VI. SOME RESULTS FROM DIFFERENT SUPERNOVA MODELS 

We have shown how the original neutrino fluxes would be modified by the MSW effect in 
a supernova even if the details of the density gradient of matter are unavailable. The input 
spectral parameters are adopted from the LL model. As compared to the LL model, the 
parameters predicted by the ongoing work of the Garching group [25,26], labeled as Gl and 
G2, suggest smaller differences in both the flux and in the average energy among neutrinos 
of different flavors: 

Gl: (E Ue ) = 12 MeV, (E De ) = 15 MeV, (E„ x ) = 18 MeV, L°jLl = L%J L° x = 0.8, 

G2: (E Ue ) = 12 MeV, (E Pe ) = 15 MeV, (E Ux ) = 15 MeV, L°jLl = L%J L° Vx = 0.5. 

To examine the impact due to the choice of models, the expected neutrino spectra for 
the six possible scenarios using the parameters of Gl and G2 are shown in Figs. 7 and 8, 
respectively. Note that the spectral shapes predicted by the LL model in Fig. 1 exhibit a 
very distinct nature from that predicted by Gl and G2. It is seen that the peak locations 
of the isopropical events for III and VI mark the most significant difference between Gl and 
G2. The peak locations for other scenarios are relatively insensitive to the choice of Gl or 
G2. As for the directional events, the six scenarios are barely distinguishable under the Gl 
model. In general, the G2 model predicts directional spectra that are sharper than that 
predicted by the Gl model. 

The inherent complexity that leads to the entanglement in the spectral shapes, the 
peak heights, and the peak locations among scenarios and models may render a thorough, 
quantitative analysis very difficult. However, as can be seen from Figs. 7 and 8, the most 
significant difference among scenarios for a given model appears in the spectra near the 
peaks. Thus, the observed i?f and R® might be able to reflect the distinct characters of 
each model. We show the plots of i?f = i?f (ip) for Gl and G2 in Fig. 9. 

It is seen that Gl gives rise to distinct values of Ri(jp) ranging from ~ 11.3 to ~ 19.9. 
As for G2, the values of R®(ip) are confined within 36 ~ 37 for ip < it/2, and 28 < i?f < 37 
for 57r/6 < ip < 7r/2. Fig. 9 suggests that if the spectral parameters of Gl are adopted, the 
six different scenarios may still cause R® to vary by roughly a factor of two. On the other 
hand for G2, the ratio i?f (VO appears to be insensitive to the uncertainties represented by 
the distinct scenarios at small incident angle. As a comparison, one recalls from Fig. 3 that 
the possible values of R®(ip) for the LL model varies from ~ 2.6 to ~ 8.8, which is in a 
region very different from that for Gl and G2. We conclude that the observable R® would 
reflect the distinct characters of the three different model. 

The ratio R$(tp) can be analyzed in a similar way. We show the plots of R^ijp) for Gl 
and G2 in Fig. 10. It is clear that the values of i?f predicted by the three different models 
fall in distinct regions: R$ = 1.2 ~ 2.5 for LL, R° = 2.7 ~ 3.5 for Gl, and Rg = 5.3 ~ 6.0 
for G2. In addition, it should be pointed out that we also perform calculations on the other 
ratios, Rf, R®, and R$, using the input parameters predicted by Gl and G2. However, 
unlike Rf and R®, the results show that the expected values for each of these ratios would 
all tangled up unless either a particular model is specified or the uncertainty due to the 
density profile can be significantly reduced. 
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In summary, the choice of models would lead to distinct signatures in the values of i?f (ip) 
and R^ip). The degeneracy among the six scenarios for either LL or Gl could be removed. 
On the other hand, the very nature of G2 would lead to only very little differences among 
the scenarios, i.e., it would be difficult for them to be quantitatively distinguished. However, 
if the correct supernova parameters are the ones predicted by the G2 model, the analysis 
based on the G2 model could be relatively free from the disturbing uncertainties, such as 
that originating from the density profile, the mass hierarchy, and the mixing angle #13. 

VII. SUMMARY AND CONCLUSIONS 

While a satisfactory knowledge of the supernova mechanism is still far from complete, 
the neutrino flux emitted from the core-collapse supernova could be a promising tool for 
probing the unknown neutrino properties. It is the intention of this work to search for pos- 
sible hints of the unknown neutrino properties using the experimental observables derived 
from the supernova neutrino signals. Given the uncertain density profile and the choice of 
input spectral models, the aim is set to exploring the unknown neutrino parameters when 
the neutrino fluxes are modified by the MSW resonant conversion in both the supernova and 
in crossing Earth. Although the MSW mechanism has been well established in the litera- 
ture, this analysis shows that the phenomenological outcomes due to the possible variation 
of density profile still provide certain new physical insight to the nature of the unknown 
neutrino properties. 

In stead of following the traditional approach based on the global density profile p ~ r~ 3 , 
or the profile with a fixed power in general: p ~ r n , this analysis allows the density profiles 
near the locations of resonance to vary independently as arbitrary power-law functions: 
p(r) = c\r ni for the lower resonance, and p(r) = c/ l r n ' 1 for the higher resonance. It is shown 
that the variation of rik {k = I or h) would cause different degrees of impact on the survival 
probability for u e (z/ e ) according to whether > —5 or < —5. As far as the survival 
probabilities are concerned, there exists six possible scenarios, which are characterized by 
possible combinations of the density profile, the mass hierarchy, and the mixing angle #13. 

We propose observable quantities that are derived from the expected supernova neutrino 
events at the water Cherenkov detector, and explore the feasibility of probing the neutrino 
parameters with these observables. As an illustration, the set of spectral parameters pre- 
dicted by the Lawrence Livermore group (LL) is adopted. To investigate how the choice of 
supernova models would affect the outcomes, the analysis based on the input parameters 
predicted by the Graching model is also presented as a comparison. 

Given the unknown density profile of the supernova and the choice of input spectral 
parameters, we show how the unknown mass hierarchy and the mixing angle #13 would leave 
observable signatures on the neutrino events. The difficulties and limits for this type of 
analysis are also discussed. It is shown that the two observables, R\ and R2, can combine 
to act as the most effective pair of observables for probing the neutrino parameters. Other 
observables, R3, R^, and R$, may be considered as the supplementary references. A more 
realistic analysis incorporating the Earth regeneration effect is also included. Figs. 3, 4, 5, 
and 6 show how each of the observables varies with the incident angle of neutrinos when the 
Earth matter induces additional flavor conversion prior to the detection of neutrinos. The 
general formulation can be further expanded to accommodate different models of the Earth 
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density profile and can be applied to analyzing the effects if the neutrinos are registered by 
other types of terrestrial detectors. 

It should be pointed out that certain non-MSW effects may occur in the complex envi- 
ronment of a supernova at different space and time scales. These effects may be independent 
of the density profile and may be totally irrelevant to the MSW conversion. If the proper 
conditions are met, recent simulations suggest that effects such as the neutrino self coupling 
in dense media (see, e.g., Ref. [34-42] and the references therein) and the neutrino flavor 
de-polarization associated with the after shock turbulence (see, e.g., Ref. [43-47] and the 
references therein) could occur. These effects are expected to impact the efficiency of neu- 
trino flavor transition. It is suggested that both the shock reheating and the hot bubble 
epochs could provide the proper conditions for two classes of large-scale collective effects: 
synchronization and bipolar flavor transformation, which arose if neutrinos themselves form 
significant background. The neutrino flavor transformation due to this neutrino coupling 
effect shares no common origin with that due to the MSW resonance. In addition, a fluctu- 
ating density in the post-shock region could be created by the turbulent convective motions 
behind the forward moving shock and the reverse shock. Such fluctuations could erase part 
of the shock wave imprint on the neutrino spectra even with relatively small amplitudes. 
It is suggested that if the amplitude of density fluctuations in the turbulence is large, the 
average neutrino survival probability would saturate to 1/2. The information about the 
initial state is then lost, and the final state is depolarized. This depolarization could leave 
unique signatures that are distinct from those left by the shock and other density effects. 

As concluded by this work, the variation of density profile, which plays a determinant role 
in the MSW flavor transition, could modify the spectra significantly and lead to variation 
of the proposed observables by a factor of two or more. Whether the flavor transition effect 
arising from one origin would dominate over that from the others depends on the spectral 
parameters, and on the neutrino intrinsic properties such as the mass hierarchy and the 
magnitude of #13. In addition, it also depends on the physical conditions that the neutrinos 
would encounter at different stages of the supernova environment and in crossing the Earth 
matter. Given the undetermined neutrino intrinsic properties and the existing astrophysical 
uncertainties, it would be intriguing to further investigate, under all possible and relevant 
physical conditions, how the non-MSW effects would act as a modification to the MSW 
effect or even dominate over the MSW effect. The more detailed analyses with reduced 
assumptions would certainly help suggest how some of the existing MSW-based paradigms 
for the neutrino flavor conversion in a supernova should be reconstructed. 
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APPENDIX A: EXPRESSIONS FOR AND ' J 

The neutrino flux of flavor a (a = e, e) for the six scenarios at the surface of Earth, F£, 
and at the detector, F^ ,J , are related to the original fluxes F® by the following: 

Fl = F e °(cos 4 6 12 + sin 4 12 ) + - (cos 4 12 + sin 4 6 12 )\, (Al) 
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F e °[cos 2 12 + P 2e (l - 2 cos 2 e 12 )\ + F°[l - (cos 2 12 + P 2e (l - 2 cos 2 2 ))], 

F/ = P°(cos 4 9 12 + sin 4 12 ) + P°[l - (cos 4 12 + sin 4 12 )], 
F°[l - (sin 2 12 + P 2e (l - 2 sin 2 12 ))] + P°[sin 2 12 + P 2e (l - 2 sin 2 12 )], 

F II _ pD,II _ pO 

Fi I = F°(cos 2 9 12 )+F°( S in 2 6 12 ), 

= P°(l - P 2e ) + F° x P 2e , 

Fl n = P e ° sin 2 9 12 + F° cos 2 12 , 

F e ^" = F e °P 2e + F°(l-P 2e ), 
pj/j _ pDjn _ F o 

± e ± e J n 

P e /y = P e ° sin 2 12 + F° cos 2 12 , 
F D,/y = F op 2e + F o (1 _p 2e); 

P/ y = F^ cos 2 0i 2 + F° sin 2 12 , 
F D,/y = F o (1 _p 2e) + F Op 2e; 

F V _ _ pO 

± e ± e - 1 x ' 

F y = P°(cos 4 12 + sin 4 12 ) + F°[l - (cos 4 12 + sin 4 12 )], 
= F^cos 2 12 - P 2e (l - 2 sin 2 12 )] + F°[sm 2 12 + P 2e (l - 2 sin 2 12 )], 

F y/ = P e °(cos 4 12 + sin 4 12 ) + F°[l - (cos 4 12 + sin 4 12 )], 
= P e °[cos 2 12 + P 2e (l - 2 cos 2 12 )] + P°[sin 2 12 - P 2e (l - 2 cos 2 12 )], 

F VI = pD,VI = pO 
e e x ' 
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APPENDIX B: FLUX VARIATION DUE TO EARTH MATTER EFFECT 



Expressions for R J e and Rg are given as the following: 

j_ sin 2 gig - cos 2 gig , 2 , , . 

1 + sin 6>i2 + cos 4 6*12 

j_ 3(sin 2 gig cos 2 gig) p 2 
3(sm i2 + cos 4 0i 2 ) + 5 



Re = 0, (B3) 



^-^^(^-- 2 M, (B4) 



■ (P 2 e-sin 2 0i 2 ), (B5) 



= 0, (B6) 



= — -^2— (P 2e -sin 2 12 ), (B7) 
l + sm t/i2 



^-^^(P 2e -sm 2 M, (B8) 



^ = (B9) 



tdV _ 3(sin 2 0i2-cos 2 0i2) = 2 
3 (sin 0i2 + cos 4 0i2 ) + 5 



1 + sm 012 + COS 4 012 



i^' = 0. (B12) 
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TABLE I. Types of scenarios that are categorized by various combinations of (P,P). Each 
of the quantity in the square bracket in the first row represents the variable that leads to the 
probability to its right. The total survival probabilities, P and P, are determined by the specific 
ni, rih, 0i3, and the mass hierarchy. 
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FIG. 1. (a) The expected isotropical spectra for the six different scenarios at a water Cherenkov 
detector, such as the Super-Kamiokande. Note that there are three groups of two-fold degeneracy: 
(I,V), (II,IV), and (III,VI). (b) The expected directional spectra for the six scenarios. Each curve 
is labeled according to the height of its peak. 
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FIG. 2. (a) The neutrinos travel through only the mantle of Earth before reaching the detector, 
(b) The neutrinos travel through the combination of mantle+core+mantle. Note that L is the total 
length of the path inside Earth and ■0 is the incident angle. 
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FIG. 3. The plot of i?f = N iso (E < 40MeV)/N iso (E > 40MeV) for all six scenarios as 
functions of the incident angle if). Note that the six scenarios form three nearly degenerate groups. 
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FIG. 4. The plot of B§ = N dir {E < 25MeV)/N dir (E > 25MeV) for all six scenarios as 
functions of the incident angle ip. 
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FIG. 5. The plot of R? = N iso (13MeV < E < 27MeV)/N dir (13MeV < E < 27MeV) for all 
six scenarios as functions of the incident angle ip. 
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FIG. 6. The plot of = (E)i SO / (E)di r for all six scenarios as functions of the incident angle if). 
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FIG. 7. The expected isotropical (a) and directional (b) spectra for the six different scenarios 
predicted by the Gl model at a water Cherenkov detector. There are three groups of two-fold 
degeneracy in (a): (I,V), (II,IV), and (III, VI). Note that the curves in (b) are barely distinguishable 
and are not labeled. 
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FIG. 8. The expected isotropical (a) and directional (b) spectra for the six different scenarios 
predicted by the G2 model at a water Cherenkov detector. There are three groups of two-fold 
degeneracy in (a): (I,V), (II, IV), and (III, VI). Each curve in (b) is labeled according to the height 
of its peak. 
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FIG. 9. The ratio J?f = iV iso (£ < 40MeF)/A7" iso (.E > 40M e y) for all the six scenarios as 
functions of the incident angle tp i n the Gl model (a), and in the G2 model (b). Note that 
11.3 < Ri < 19.9 in (a), and that the six scenarios are almost indistinguishable for tp < ir/2 in 
(b), 36 < R? < 37. 
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FIG. 10. The ratio B§ = N dir (E < 25MeV)/N dir (E > 25MeV) for all the six scenarios 
as functions of the incident angle ip in the Gl model (a), and in the G2 model (b). Note that 
2.65 < J?f < 3.50 in (a), and that R? > 5.3 in (b). 
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